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Abstract 

The equilibrium of a resistive axisymmetric plasma with purely toroidal 
flow surrounded by a conductor is investigated within the framework of the 
nonlinear magnetohydrodynamic theory. It is proved that a) the poloidal 
current density vanishes and b) apart from an idealized case the pressure 
profile should vanish on the plasma boundary. For the cases of isothermal 
magnetic surfaces, isentropic magnetic surfaces and magnetic surfaces with 
constant density, the equilibrium states obey to an elliptic partial differ- 
ential equation for the poloidal magnetic flux function, which is identical 
in form to the corresponding equation governing ideal equilibria. The con- 
ductivity, which can not be either uniform or a surface quantity, results 
however in a restriction of the possible classes of equilibrium solutions, e.g. 
for the cases considered the only possible equilibria with Spitzer conduc- 
tivity are of cylindrical shape. 
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I. Introduction 



The understanding of the equilibrium properties of a magnetically confined 
plasma is one of the basic objectives of fusion research. The major part of the 
studies to date has been devoted to static magnetohydrodynamic (MHD) equi- 
libria of an axisymmetric plasma. They are governed by a partial elliptic dif- 
ferential equation for the poloidal magnetic flux function ifj, which contains two 
surface quantities [i.e. the pressure P(ip) and the poloidal current I{ip)], known as 
Grad-Schliiter-Shafranov equation (Lust and Schluter 1957; Grad 1958; Shafra- 
nov 1958). 

In the presence of flow the equilibrium becomes more complicated. Although 
equilibrium studies of flowing plasmas began in the mid 1950s (e.g. Grad 1960), 
since the early 1980s there has been an increasing interest in the investigation of 
the equilibrium properties of plasmas with mass flow (Zehrfeld and Green 1972; 
Maschke and Perrin 1980; Morozov and Solovev 1980; Hameiri 1983; Clemente 
and Farengo 1984; Semenzato et al 1984; Kerner and Tokuda 1987; Throumoulopou- 
los and Pantis 1989; Avinash et al 1992; Zelazny et al 1993; Throumoulopoulos 
and Pantis 1996; Throumoulopoulos and Tasso 1997) which was motivated by the 
observation of plasma rotation in many tokamaks heated by neutral beams (Suck- 
ewer et al 1979; Brau et al 1983; Core et al 1987; Scott et al 1989; Tammen et al 
1994) With the adoption of a specific equation of state, e.g. isentropic magnetic 
surfaces (Morozov and Solovev 1980), the axisymmetric ideal MHD equilibrium 
states obey to a partial differential equation for tp, containing five surface quan- 
tities, in conjunction with a nonlinear algebraic Bernouli equation. Unlike the 
case in static equilibria, the above mentioned differential equation is not always 
elliptic; there are three critical values of the poloidal flow at which the type of 
this equation changes, i.e. it becomes alternatively elliptic and hyperbolic. To 
solve the equilibrium problem in the two elliptic regions, several computer codes 
have been developed (Semenzato et al 1984; Kerner and Tokuda 1987; Zelazny et 
al 1993). There is experimental evidence, however, that either the flow is purely 
toroidal (Suckewer et al 1979; Core et al 1987; Scott et al 1989; Tammen et al 
1994) or its poloidal component is one order of magnitude lower than the toroidal 
one (Brau et al 1983). It was also found that the poloidal flow is efficiently damped 
by magnetic pumping (Hassam and Kulsrud 1978). In addition, toroidal rotation 
has been suggested as a possible means of creating a magnetic field in the sun 
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and stars (Plumpton and Ferraro 1955) and of stabilizing the tilting (Hayashi 
and Sato 1987), ballooning (Miller et al 1995) and drift-like (Sen and Rusbridge 
1995) modes. Therefore toroidal flow may be considered a dominant effect. For 
a purely toroidal velocity component the partial differential equation becomes 
elliptic and can be solved analytically (Maschke and Perrin 1980; Clemente and 
Farengo 1984; Throumoulopoulos and Pantis 1989). 

The investigation of the general MHD equilibrium, involving plasma flow, 
finite conductivity and viscosity is a very difficult problem, particularly in light 
of the uncertainties in the viscous stress tensor. For example, for the case in 
which the mean free path is greater than a gyroradius, as it is in the case of a 
tokamak, the "ion parallel viscosity" (Braginskii 1965; Balescu 1988) takes an 
extraordinary large value, so large as to make the applicability of either version 
of the MHD viscosity stress tensor to tokamak dynamics questionable. It may 
be noted, however, that the ion parallel viscosity first increases strongly with 
the temperature until the mean free path reaches the order of magnitude of the 
machine, then decreases for larger values of the mean free path. In particular, 
in the high viscosity regime resistive equilibria with flows originated from paired 
toroidal vortices, which are very likely created in a toroidal magnetofluid, were 
investigated recently (Montgomery Bates and Li 1997a); the existence of such 
vortices is however conjectured to characterize magnetofluids beyond the high- 
viscosity limit. Thus, the form and the magnitude of the appropriate viscous 
stress tensor to be used in tokamak MHD has been the subject of considerable 
discussion, and at present these discussions show no signs of converging. 

For an inviscid plasma Pfirsch and Schliiter (Pfirsch and Schliiter 1962; see 
also Wesson 1987) showed long time ago that in the collisional regime the toroidal 
curvature gives rise to an enhanced diffusion, which is related to the conductivity 
parallel to the magnetic field. Since then there have been attempts to examine 
the effects of finite conductivity (e.g. Grad and Hogan 1970), but these attempts 
appear not having led to a satisfactory resolution; too many possibilities are 
raised to be able to deal with any one of them conclusively. In the above men- 
tioned studies the inertial-force flow term p(v • V)v is neglected in the equation 
of momentum conservation. For ion flow velocities of the order of 10 7 cm/sec, 
which have been observed in neutral-beam-heating experiments (Suckewer et al 
1979; Brau 1983 et al; Tammen et al 1994), however, the term p(v • V)v can not 
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be considered negligible. Therefore, it is worthwhile to investigate the nonlinear 
resistive equilibrium. Such an investigation is also encouraged by a proof (Tasso 
1979) according to which a conductivity profile remaining constant on a magnetic 
surface is not compatible with the Grad-Schliiter-Shafranov equation. The non 
existence of static axisymmetric resistive equilibria with a uniform conductivity 
was also suggested recently (Montgomery and Shan 1994; Bates and Lewis 1996; 
Montgomery et al 1997b). 

In the present paper we study the MHD equilibrium of an axisymmetric 
plasma with scalar conductivity and purely toroidal flow surrounded by a con- 
ductor, including the term p(v • V)v in the momentum conservation equation. 
It is found that the presence of flow does not remove the above mentioned in- 
compatibilities associated with the conductivity. Also, the conductivity has an 
impact on the direction of the current density, on the boundary conditions and 
restricts the possible classes of equilibrium solutions. 

The paper is organized as follows. For the system under consideration equilib- 
rium equations independent of the plasma equation of state are derived in Sec. II. 
The cases of equilibria with (a) isothermal (b) isentropic and (c) constant density 
magnetic surfaces are examined in Sec. III. Sec. IV summarizes our conclusions. 

II. Equilibrium equations 

The MHD equilibrium states of an inviscid plasma with scalar conductivity 
are governed by the following set of equations, written in standard notations and 
convenient units: 

V • (pv) = (1) 
p(v V)v=j x B- VP (2) 

V x E = (3) 
VxB=j (4) 

V • B = (5) 

E + vxB = l (6) 

a 

It is pointed out that, unlike to the usual procedure followed in equilibrium stud- 
ies with flow (Zehrfeld and Green 1972; Morozov and Solovev 1980; Hameiri 
1983; Semenzato Gruber and Zehrfeld 1984; Kerner and Tokuda 1987; Zelazny 
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et al 1993) in the present work an equation of state is not included in the above 
set of equations from the outset and therefore the equation of state independent 
Eqs. ( |T7D and (|T%| ) below are first derived. This alternative procedure is conve- 
nient because the equilibrium problem is further reduced then for specific cases 
associated with several equations of state. 

The system under consideration is a toroidal axisymmetric magnetically con- 
fined plasma, which is surrounded by a conductor (Fig. 1). Cylindrical coordi- 
nates R, 0, z are employed with unit basis vectors e#, e^, e z , where e z is parallel 
to the axis of symmetry. The position of the surface of the conductor is speci- 
fied by some boundary curve in the (R, z) plane. The equilibrium quantities do 
not depend on the azimuthal coordinate <fi. Also, it is assumed that the plasma 
elements flow solely along e^. Eq. (HI) is then identically satisfied. 

For the above mentioned system we first prove that the equilibrium poloidal 
current density must vanish. Inserting the expression 

B = + B pol = B^ + ^Vip x (7) 
for the magnetic field into Ampere's law @) yields for the current density: 

j = --^A*V>e* + -^V(iLB^) x e^. (8) 

Here, tp is the poloidal magnetic flux function which labels the magnetic surfaces, 
and A* is the elliptic operator defined by Aijj = i? 2 V ■ (Vtp/R 2 ). With the aid 
of Eqs. (0) and (§), the component of the force balance equation @ along is 
written in the form 

e ■ V(i^) x W> = 0. (9) 

Eq. (H) implies that RB^ is a surface quantity: 

RBf = I(ip). (10) 

Integration of Ohm's law (||) along a contour c defined by the cut of an arbitrary 
closed magnetic surface with the poloidal plane (Fig. 1) leads to the equation 

f-E-d\+ /(vxB).dl= l-{yd\). (11) 

Since in equilibrium it holds that dH/dt = 0, the first integral vanishes by Stokes's 
theorem. Also, the second integral vanishes because v is purely toroidal. Con- 
sequently, the term on the right-hand side must vanish, implying that jp Q j = 
and therefore I = constant. 



6 



The 0-component of the electric field associated with the current density can 
be obtained from Faraday's law (^]): 

E* = & = £7 §, (12) 

where E is the value of the electric filed in the cylindrical surface R = Rq. An 
electric field of this kind can be created by a linear change in time of the magnetic 
flux confined to an axisymmetric very long iron core, which is directed along e z 
and is located within the hole of the torus. This can be considered as a highly 
simplified representation of a tokamak. 



Eq. (J12 ) has an impact on the boundary conditions, i.e. the component of E 
tangential to the plasma-conductor interface does not vanish. Therefore, apart 
from a special idealized situation in which the surface of a perfectly conducting 
container is coated by a very thin layer of insulating dielectric (Bates and Lewis 
1996; Montgomery et al 1997b) the container can not be considered perfectly 
conducting. Accordingly, Ohm's law with finite conductivity applied in the vicin- 
ity of the plasma-conductor interface does not permit the existence of a surface 
layer of current. It is now assumed that the position of the conductor is such 
that its surface coincides with the outermost of the closed magnetic surfaces. 
Thus, the condition B • n = 0, where n is an outward unit vector normal to the 
plasma surface (Fig. 1), holds in the plasma-conductor interface. It is noticed 
that this is possible only in equilibrium, because in the context of resistive MHD 
time dependent equations, the magnetic flux is not conserved. Integration of 
the momentum conservation equation across the plasma-conductor interface and 
projection of the resulting equation parallel to n yields (with B being continuous 
in the plasma-conductor interface) 

P + p(n-v) 2 = 0. (13) 

Since all the terms of Eq. (0) are non-negative and the flow is purely toroidal, 
it follows that the pressure P must vanish at the plasma boundary. 

To reduce the equilibrium equations ([j])-(|6]), a second surface quantity is iden- 
tified by projecting the equation 

V x (v x B - = (14) 
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[following from a combination of Eqs. (|3|) and (P)] along e^,. This yields 



(15) 



Equation (p~5|) implies that the angular frequency of the flow shares the same 
surfaces with ip: 

W = ^ = u;(V0. (16) 

With the aid of equations (0), (|8|) with / = const., flTU|) and flTB|), the components 
of Eq. (0) along Bp Q j and perpendicular to a magnetic surface are put in the 
forms, respectively, 



and 



B pol 

A*<0 + # 







(17) 



2 vv> 



VP - pcu 2 V 



i? 2 



0. 



It is noted here that Eqs. (|9]), ( |i~7|) and (|T8| ) are orthogonal components of Eq. 
(|2D and therefore they are independent. 

Summarizing, the resistive MHD equilibrium of an axisymmetric toroidal 
plasma with purely toroidal flow is governed by the set of Eqs. (fi7|), [18]) and 



A*i/j + RqE q o = 0, 



(19) 



following from Eqs. (§) and (|I 
flow, which leads to jp Q j = and uj 



Owing to the purely toroidal direction of the 



18) do not contain 



17D and 

the conductivity and are identical in form to the corresponding equations govern- 
ing ideal equilibria. Therefore, on the one side, several properties of the ideal 
equilibria, e.g. the Shafranov shift of the magnetic surfaces and the detachment 
of the isobaric surfaces from the magnetic surfaces (see the discussion following 
Eq. (p2|) below) remain valid. On the other side, Eq. (O) can be employed as 



an algebraic equation for the conductivity i.e. the solutions of the "ideal MHD" 
equations ([[7]) and (|l^) can determine the spatial dependence of a by means of 
the formula 

' -Aif>*. (20) 



a 



RqEq 



Eq. (|20| ) implies that the conductivity can not be uniform, as the case in static 
equilibria (Tasso 1979; Montgomery and Shan 1994; Bates and Lewis 1996; Mont- 
gomery et al 1997b). 
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It may be noted here that the equilibrium problem can be treated in a different 
way, suggested for the case of static resistive equilibria (Bates and Lewis 1996), 
i.e. Eq. (fKj| ) is employed as an elliptic partial differential equation for ip, which is 
solved after the spatial dependence of a is assigned. To our opinion, this method 
has two drawbacks: a) A physically unjustifiable assignment of conductivity is 
required from the outset, b) The conductivity profile must be compatible with Eq. 
(|i"8|), a requirement which (at least in the presence of flow) makes the assignment 
difficult. This becomes more clear in next section [see Eqs. (|23|), (|29|) and ([331)1. 

Ill Magnetic surfaces with special properties 

In order to further reduce the equilibrium equations (|lj]) and (|Tj^) and find 
out their influence on the conductivity profile, the starting set of equations (|l])-(||) 
must be supplemented by an equation of state, e.g. P = P(p, T), along with an 
equation determining the transport of internal energy such as 

° v (J^ + V ' V ) T + PV - v + kV2t = > ( 21 ) 

the energy transport equation for a perfect gas, where C v is the volume specific 
heat, k the thermal conductivity and = 0. Such a rigorous treatment, how- 
ever, makes the equilibrium problem very cumbersome. Alternatively, since for 
fusion plasmas the thermal conduction along B is expected to be fast in relation to 
the heat transport perpendicular to a magnetic surface, equilibria with isothermal 
magnetic surfaces is a reasonable approximation (e.g. Mashke and Perrin 1980; 
Clemente and Farengo 1984; Throumoulopoulos and Pantis 1989). In particular, 
the even simpler case of isothermal resistive equilibria (T = constant) has also 
been considered (Grad and Hogan 1970). Equilibria with isothermal magnetic 
surfaces are examined as follows. 

Equilibria with isothermal magnetic surfaces 

In addition to T = T(ip) the plasma is assumed to obey the perfect gas law 
P = XpT. For this kind of equilibria, Eq. (|D]) can be integrated yielding an 
expression for the pressure: 

P = F r (i/,)exp(^\. (22) 
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We note here that, unlike in static equilibria, in the presence of flow magnetic 
surfaces do not coincide with isobaric surfaces because Eq. (H) implies that B- VP 
in general differs from zero. In this respect, the term Pt{iP) is the static part of 
the pressure (u 



0). Inserting Eq. (22) into Eq. 



A*ijj + R 2 



Pt + Pi 



R 2 




exp 



T8T), the latter takes the form 
R 2 u 2 ^ 



2XT 



0. 



(23) 



where the prime denotes differentiation with respect to ip. This is the equilibrium 
equation for a resistive plasma with purely toroidal flow and isothermal magnetic 
surfaces. It contains the three arbitrary surface quantities T(ip), u)(if>) and Pr(if>), 
which must be found from other physical considerations. For u = it reduces to 
the Grad-Schliiter-Shafranov equation for vanishing poloidal current. For ideal 



equilibria, Eq. (23) was obtained by Maschke and Perrin (1980). For special 
choices of the surface quantities, which must lead to physically reasonable profiles 
for the equilibrium quantities, Eq. (^) can be linearized and solved analytically. 
In particular, with the ansatz u 2 /(2XT) = lot = const, it takes the simpler form 



A*V + P 2 PrWexp (R} 



IjJT 



(24) 



Solutions of Eq. (p4[) have been obtained for P? oc if> (Maschke and Perrin 
1980) and for P? oc if) 2 (Clemente Farengo 1984; Throumoulopoulos and Pantis 
1989) and the detachment of the isobaric surfaces from the magnetic surfaces was 
evaluated. 

Eqs. (0) and (H) yield 



EqRqu = R 



P T + P: 



R 2 




■R 2 uj 2 " 
2XT 



exp 



(25) 



As the case in static equilibria (Tasso 1970; Montgomery and Shan 1994) toroidal 
equilibria with the conductivity remaining constant in a magnetic surface are not 
possible. Indeed, if a = a(ip) the left-hand side of Eq. (|25|) becomes a surface 
quantity, while the right-hand side (in addition to surface quantities) involves an 
explicit dependence of R. Considering the variables R and ip, instead of R and 
z, as independent [then z = z(R,i[))], Eq. (p5|) , after expanding the exponential 
term in a Taylor series, can be written in the form J2^=o Uni^R™ — 0. The 
latter requires that the relation a n (ip) = must be satisfied for all n. This is 
impossible because a = E R a(ip) and therefore it follows that if) = if>(R). Thus, 
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the only possible a = <j(ip) equilibria are of cylindrical shape. It may be noted 
here that the non existence of another class of resistive toroidal MHD equilibria, 
i.e. axisymmetric incompressible f3 p — 1 equilibria with purely poloidal flow, was 
proved by Tasso (1970). 

To find out whether the above incompatibilities associated with the conduc- 
tivity profile arise from the choice T = T(ip), one alternatively can consider the 
cases of isentropic processes and incompressible flows, which has been the sub- 
ject of extensive ideal MHD studies and are, respectively, associated with isen- 
tropic magnetic surfaces (Maschke and Perrin 1980; Morozov and Solovev 1980; 
Hameiri 1983; Semenzato et al 1984; Kerner and Tokuda 1987; Throumoulopou- 
los and Pantis 1989; Zelazny et al 1993) and constant density magnetic surfaces 
(Avinash et al 1992; Throumoulopoulos and Pantis 1996; Throumoulopoulos and 
Tasso 1997). 

Equilibria with isentropic magnetic surfaces 

We consider a plasma with large but finite conductivity such that for times 
short compared with the diffusion time scale, the dissipative term m j 2 /cr can 
be neglected. This permits one to assume that the magnetic surfaces are isen- 
tropic, i.e S = S(ip) where S is the plasma specific entropy. It may be noted 
that S = S(ip) can be proved for the stationary states of an ideal plasma with 
arbitrary flows on the basis of the entropy conservation (Morozov and Solovev 
1980; Hameiri 1983). This case was also studied for ideal equilibria with purely 
toroidal flow (Maschke and Perrin 1980; Throumoulopoulos and Pantis 1989). 
In addition, the plasma is assumed to being a perfect gas whose internal energy 
density E is simply proportional to the temperature. Then, the equations for the 
thermodynamic potentials lead to (Maschke and Perrin 1980) 

p = A{S)p 1 (26) 

and 

E= m^ = a. (27) 

7- 1 7 

Here, A = A(S) is an arbitrary function of the specific entropy S = S(tp), 
H = E + P/p is the specific enthalpy and 7 is the ratio of specific heats. For 
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simplicity and without loss of generality we choose the function A to be identical 

(|17D yields 

H- — = 9(iP), (28) 



with S. Consequently, integration of Eq. (|T7D yields 

R 2 lu 2 



where 9 is a surface quantity. Eq. (|T8| ) reduces then to 
A^ + R 2 (l + 



2tf 



x 



5 



i? 2 



s 



1-7) 



LULU 



'6 s 



( S l- V y 



to 
26 



0, (29) 



where r] = 7/(7 — 1). This is the resistive equilibrium equation for isentropic 
magnetic surfaces, which is the same as the ideal equation (Maschke and Perrin 
1980; Throumoulopoulos and Pantis 1989) with jp Q j = 0. The three surface 
quantities are now u)(if)), S(ip) an d 9{ip). With the ansatz lu 2 /(29) = lu s = const., 
Eq. (|29|) is put in the simpler form 



A^ + R 2 (l + u s R 2 ) V {p s y = 0, 



(30) 



where the function 



'9" 



S 



i-n 



stands for the static equilibrium pressure. Solutions of Eq. ( p0|) are avail- 
able in the literature for P s oc ip (Maschke and Perrin 1980) and for P s oc ip 2 
(Throumoulopoulos and Pantis 1989). 

Eqs. (|20|) and (|29|) imply that toroidal isentropic equilibria are not compatible 
with er = a(ip). 111 addition, equilibria with Spitzer conductivity a = aT 3 ^ 2 satisfy 
the equation 



RqEqCL 



X 



1 

7^ 



R 2 lu t 



9 + 



3/2 



S 



1-7) 



+ R 2 



-- R 2 \1 + 



R 2 lu 2 V- 1 



S 



1 — 7) 



29 

lulu' 



9' 



+ 



( S l- V y 



to 
29 



0,(31) 



following from Eqs. (|20| ) and (p6|-^9|). Eq. (|31~D contains, in addition to surface 
quantities, the coordinate R explicitly. Therefore, as in the case of isothermal 
magnetic surfaces and a = <j(ip), equilibria with isentropic magnetic surfaces and 
Spitzer conductivity must have cylindrical shape. 
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Equilibria with magnetic surfaces of constant density 



Equilibria with incompressible flows are usually associated with constant den- 
sity magnetic surfaces (Avinash et al 1992; Throumoulopoulos and Pantis 1996; 
Throumoulopoulos and Tasso 1997). For axisymmetric equilibria with purely 
toroidal flow, Eqs. (p]) and V • v = are satisfied identically and provide no 
information about p. Nevertheless, p = p(ip) is a physically possible equation of 
state. For this case, integration of Eq. flTTj) leads to 



pR 2 co 2 



P = P P W) + 



(32) 



where P p (if)) is the static equilibrium pressure. Consequently, Eq. ( 18|) reduces 
to 

I? 
2 



A*<0 + R 2 



0. 



(33) 



As in the cases of isothermal and isentropic magnetic surfaces, Eqs. ( p0|) and (|33|) 
imply that in toroidal equilibria the conductivity can not be a surface quantity. 
Also, equilibria with Spitzer conductivity satisfy the relation 



E R a 



Xp 



+ 



u 2 R 2 \ 3/2 



2A 



R 2 



P> + * 
p 2 



puj' 



(34) 



which implies that magnetic surfaces must have cylindrical shape. 

Concluding this subsection, we derive two classes of exact solutions of Eq. 
(|33|) . With the use of the ansatz 



pu) 



const, 



(35) 



Eq. (|33| ) reduces to 



A*V> + R 2 P' P = 



(36) 



which is identical in form to the the Grad-Schluter-Shafranov equation for van- 
ishing poloidal current. A simple solution of Eq. (|36|) can be derived for P p oc ip 
(Shafranov 1966): 



i> = ^ c % (2R 2 c -R 2 -4d 2 z 2 



Rt 



(37) 



Here, ip c is the value of the flux function at the position of the magnetic axis (z = 
0, R = R c ) and d 2 is a parameter which determines the shape of the flux surfaces, 
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e.g. for d 2 = 1 Eq. fl3"T|) describes a compact torus contained within a spherical 
conductor. For P p oc -?/> 2 , solutions of Eq. ([3(J) separable in the coordinates R and 
z can be expressed in terms of Coulomb wave functions and describe equilibria 
of either a tokamak with rectangular cross section (Hernegger 1972; Maschke 
1973) or a compact torus contained within a rectangular conductor (Berk at al 
1981). These solutions can yield a peaked current density profile vanishing on 
the boundary. Owing to the presence of the flow term in Eq. (|32|) however the 
pressure does not vanish on the boundary, as it should for resistive equilibria. 
Thus, solutions of this kind are not valid in the vicinity of the boundary. The 
drawbacks of the solutions presented above in this paragraph, arise from the 
choices (|35|) and P p oc ip. They can be eliminated if, alternatively, the ansatz 
puj 2 = kip 2 and P p = fiip 2 /2 is employed, where k and fi are constants. Then, Eq. 
( |33|) reduces to 

A*V> + R 2 [n + kR 2 ] if> = 0. (38) 

The last equation can not be solved analytically and numerical solutions will not 
be derived here. 

IV. Conclusions 

The equilibrium of an axisymmetric plasma with scalar conductivity and 
purely toroidal flow surrounded by a conductor was investigated within the frame- 
work of the nonlinear MHD theory. It was found that a) the poloidal current 
density jp Q | must vanish and b) apart from an idealized case, perfect-conducting- 
boundary conditions are not relevant. The equilibrium states are governed by a 
set of two differential equations [Eqs. ( |TTD and fllS])], which do not contain the 
conductivity and are identical in form to the corresponding equations for ideal 
stationary equilibria. The conductivity, however, restricts the allowed equilibria 
to those that satisfy jp Q j = and have pressure profiles vanishing on the bound- 
ary. With the assumption of a specific equation of state, the above mentioned 
set of equations reduces to a single elliptic partial differential equation for the 
poloidal magnetic flux function which contains three surface quantities [e.g. Eqs. 
(|23"D, (|29"D , (|33"D 1. The cases of isothermal magnetic surfaces, isentropic magnetic 
surfaces and magnetic surfaces of constant density were examined and, for spe- 
cific choices of the surface quantities involved in each case, equilibrium solutions 
were presented and discussed. 
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As the case in static equilibria, the conductivity can not be uniform. In addi- 
tion, for all three cases of equation of state examined, the conductivity can not 
be a surface quantity; the only possible equilibria with conductivity remaining 
constant in a magnetic surface have cylindrical shape. Also, for isentropic mag- 
netic surfaces and magnetic surfaces of constant density, it was shown the non 
existence of toroidal equilibria with Spitzer conductivity. 

It is interesting to investigate whether these unusual properties of conductivity 
persist when additional physical input, e.g. a general flow with non vanishing 
poloidal and toroidal components is included. Of course, pursuing such kind 
of investigations, one gets a step closer to the cumbersome problem of MHD 
equilibrium with flow, finite conductivity and viscosity. The results of the few 
studies which have been conducted to date indicate, however, that this objective 
may be of crucial importance for fusion plasmas. 
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Figure 



FIG 1 The geometry of an axisymmetric toroidal plasma surrounded by a con- 
ductor. 
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